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ABSTRACT 

This is an introduction to the equation of state (EOS) for stellar matter at subnuclear densities 
calculated with the Statistical Model for Supernova Matter (SMSM). Tables for EOS of the supernova 
matter are constructed for a wide range of the densities p = (10~® — 0.32)/9o, and of the temperatures 
T = 0.2 — 25 MeV for different electron fractions Ye = 0.02 — 0.56. The cases of /3-equilibrium are also 
considered. In this paper we explain the data structure and content of the SMSM EOS tables. 
Subject headings: equation of state - nuclear astrophysics, abundances - stars: neutron - supernovae 



1. INTRODUCTION 

As established in intensive experimental studies of last 
20 years, many nuclear reactions lead to the formation of 
thermalized nuclear systems characterized by subnuclear 
densities and temperatures of 3-8 MeV. De-excitation 
of such systems goes through nuclear multifragmen- 
tation, i.e. break-up into many excited fragments 
and nucleons. As it is generally accepted, thermal 
and chemical equilibrium can be established in many 
multifragmentation reactions. Transport theoretical 
calculations of both central heavy-ion collisions around 
the Fermi-energy and peripheral heavy-ion collisions 
predicts momentum distributions of nucleons which are 
similar to equilibrium ones after f irst ~ 100 fm/c (see 
references in the topical issue of lEur. Phvs. J. A. 30 I 
(|2006j ). The final proof of equilibration was provided 
by numerous comparisons of statistical models with 
experimental data: These models perfectly describe 
many characteristics of nuclear fragments observed 
in multifragmentation experiments: multiplicities of 
intermediate-mass fragments, charge and isotope dis- 
tributions, event by event correlations of fragments 
(including fragments of different sizes), their angular 
and v el ocity correlation s , and o ther observables (iGross I 
p9^; 'Bondorf et al. " (1995^; 'Souliotis et al. " (20071; 
lOgul et al. (2011); Botvina, Iliinov. & Mishustin 

(|l990|) : iBotvina et al. I (119951) : Scharen berg et al. 

(12001): iD^gostino et all (119961 jTgQa): .BeUaize et al 



200l: lA^evev et al. I (I2002D: iHauger et al. I (|2000|): 



llglio et al. ' C2006'): 'Hudan et al. I (I2009D: iWarig et al. I 
Jl999); Viola ctal. (20011): iRodionov et aL I (|2002t ): 
[Pienkowski et al. I ()2002[ )). The parameters of nuclear 
matter at the moment of formation of hot fragments, 
such as temperature and density, can also be estab- 
lished reliably in experiment by using observed relative 
velocities of fragments, and relative isotope yields. 

A nuclear matter at similar conditions is expected in 
many astrophysical situations, in particular, in super- 
novae. In order to clarify the thermodynamical condi- 
tions obtained in nuclear reactions, we show the phase 
diagram for symmetric and asymmetric nuclear matter 



in Fig. [T] (jBuvukcizmeci et al. I ()2013[ )). for a range of 
densities and temperatures appropriate for core-collapse 
supernova explosions. Typical conditions associated with 
multifragmentation reactions are indicated by the shaded 
area in Fig. [T] These reactions give us a chance to access 
hot nuclei in the environment of other nuclear species in 
thermodynamical equilibrium as we expect in supernova 
matter. The properties of these nuclei can be directly 
extracted from experimental data and then this informa- 
tion can be used for more realistic calculations of nuclear 
composition in stellar matter. This is a new possibility to 
study this matter besides theoretical approaches which 
use nuclear forces extracted from experimental study of 
nuclear structure. As one can estimate from Fig.[l] in the 
course of massive star collapse the stellar nuclear mat- 
ter passes exactly through the multifragmentation region 
with typical entropy per baryon S/ B = 1 — 4. The elec- 
tron fraction Yg, which is equal to the total proton frac- 
tion, in the supernova core varies from 0.1 to 0.5. In 
the two-phase coexistence region (below the solid and 
dot-dashed lines) at densities p « 0.3 — O.Spo the matter 
should be in a mixed phase, which is strongly inhomo- 
geneous with intermittent dense and dilute regions. In 
the coexistence region at lower densities, p < 0.3po the 
nuclear matter breaks up into compact nuclear droplets 
surrounded by nucleons. These relatively low densities 
dominate during the main stages of stellar collapse and 
explosion. Under such conditions one can expect a mix- 
ture of nucleons and light and heavy nuclei. In Fig. [T]we 
demonstrate also isentropic trajectories in nuclear mat- 
ter (dashed curves) and the trajectories of density and 
temperature inside the supernova core, which are taken 
from th e supernova simula t ion o f the star with 15 solar 
masses (jSumivoshi et al. I ([20,05!)) (red dotted curves). 
The snapshots in the supernova dynamics are selected 
for three stages before, at, and after the core bounce. 
In the gravitational collapse from the iron core of the 
star (BB: just before bounce), the density and tempera- 
ture roughly follows the isentropic curves S/B=l-2. At 
the core bounce, when the central density increases just 
above the nuclear matter density pQ, the temperature of 
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with the hydro-dynamical instabihties in 2D and 3D 
supernova simulations. 

2. STATISTITAL MODEL FOR SUPERNOVA 
MATTER (SMSM) 

The Statistical Model fo r Supernova Matter (SMSM 
was fi rst introduced in Ref. (jBotvina fc Mishustin I (|200 
I20l 0f)) as generalization of the Statistical Multifragmen- 
tation Model (Bondorf et al. (1995)) which was success- 
fully used for description of nu c lear multifragm e ntatiq n 
reactions (lOgul et al I (120111): iBotvina et al. I (1199 " 
Scharcnbe rg et al. I (I2001D: iD'Agostino et aLl HM' 



1999i): iBell aize et al. ' ^2002^; 'Avdevev et al. ' (2002 
Hauger et a l. (2000); Iglio et al. (2006); Hudan et al." 
( 20091) : IWang et al. I (I1999D : iViola et al. I (120011) '). This 
gives us a confidence that this model can realize the ap- 
proach to diluted stellar matter described in Introduc- 
tion. We describe supernova matter as a mixture of nu- 
clear species, electrons, photons, and perhaps neutrinos 
in thermal equilibrium. Below we give the description of 
model with latest modifications. 



Fig. 1. — Nuclear phase diagram in the 'temperature - 
baryon density' plane. Solid black and dashed-dotted pur- 
ple lines indicate bou ndaries of th e liq uid-gas coexistence 
region for symmetric ()Lamb et al. I ((1978)) and asymmet- 
ric matter (Ye = 0.2) calculated with TMl interactions 
(|Sugahara fc TokD (|T99j)). The shaded area corresponds to 
typical conditions for nuclear multifragmentation reactions. 
The dashed black lines are isentropic trajectories character- 
ized by constant ent ropy per baryon, S/B =1, 2, 4 and 6 cal- 
culated with SMSM (|Botvina fc Mish ustin 1 (120101') 1. The dot - 
ted red lines show model results of Sumiyo shi et al. 1 (|2005l ) 
for BB (just before the bounce), CB (at the core bounce) and 
PB (the post bou nce) in a core-collapse supernova. This fig- 
ure is taken from IBuyukcizmeci et al. I (|2013l ) (Color version 
online.). 



the inner core becomes higher than 10 MeV due to the 
passage of the shock wave (CB: core bounce). The tem- 
perature of the whole supernova core is still high at 150 
ms after the core bounce (PB: post bounce). The shock 
wave is stalled around 130 km in this ID calculation, 
which does not lead to an explosion. 

These conditions inside the supernova core pass 
through interesting regions of the phase diagram. The 
BB and CB trajectories go through the multifragmenta- 
tion region, which motivates us to assess the supernova 
conditions by the multifragmentation reactions taking 
place in heavy-ion collisions of intermediate energy. The 
trajectories after the bounce (CB and PB) traverse the 
phase boundary between the mixture of nuclei and the 
nucleon gas. This regio n is dominated by ligh t nucle i 
C^He and lighter ones) (IBotvina &: Mishustin I (120041): 



20101): iHempel et al.1 
jHcmpel et al. I (1201 IT 
|Ar 



LArco: 



IHempel fc Schaffner-Biehch I 
(I2012D: iTvpel et al. I mii 
iSumivoshi fc Ropke I (|2008D 
Since the dynamics of the shock wave is affected by the 
interaction of neutrinos with nucleons and nuclei, it is 
important to determine the composition of hot dense 
matter in this region by involving into consideration the 
full en semble of nuclear species ([Botvina fc Mishustin I 
(|2004f )). In fact, the heating through neutrino absorp- 
tion contributes to the revival of shock wave together 



2.1. Equilibrium conditions 

We assume that each particle i with baryon number 
Bi, charge Qi and lepton number Li is characterized by 
a chemical potential /i^, which is represented as 



fli = B.fiB + QifJ-Q + LifiL- 



(1) 



Here /is , ^q and fj.^ are three independent chemical po- 
tentials which are determined from the conservation of 
total baryon number B = J^i Bi electric charge Q — 
y^^ Qi and lepton number L = ^^ Li of the system. The 
chemical potentials for nuclear species {A,Z), electrons 
(e^je"*"), and neutrinos (i/, ?■) are expressed as 



I-IAZ = AfJ-B + ZflQ , 



(2) 



These relations are valid also for nucleons, /i„ = ^b and 
Mp — M-B + MQ • The corresponding conservation laws can 
be written as 



B 



P 



F-E^" 



V 



AZ 



AZ 



PQ 



Q 



= 77 =Y1 ^P^^ - Pe = 



^ AZ 

PL = Pe+ Pu - Pv 



(3) 



YlP 



Here p^ = p^- — p^+ is the net electron density, Y^ is the 
lepton fraction, and V is the normalization volume. The 
second equation requires that any macroscopic volume of 
the star (V) is electrically neutral. 

The lepton number conservation is a valid concept only 
if V and v are trapped in the s ystem within the neutrino 
sphere (P rakash et al. I ()1997D ). If they escape freely 
from the star, the lepton number conservation is irrel- 
evant and PL = 0. In this case two remaining chemical 
potentials are determined from the conditions of baryon 
number conservation and electro-neutrality. Since the 
/3-equilibrium may not be achieved in a fast explosive 
process, we also often fix the electron fraction Y^, in the 
calculations. 
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2.2. Ensemble of nuclear species 

The Grand Canonical v er sion of the SMM 
(jBotvina. Iliinov. fc Mishustin I (|1985t )l is used for 
describing an ensemble of nuclear species under super- 
nova conditions. After integrating out translational 
degrees of freedom one can write the pressure of nuclear 
species as 



-Pnuc = Ty^pAZ 



(4) 
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^7^ {Paz - i^Az) 



where pAz is the density of nuclear species with mass A 
and charge Z. Here gAZ is the ground-state degeneracy 

1 /9 

factor of species {A, Z), At = {2-Kh? /mfqT^ is the nu- 
cleon thermal wavelength, rajsi « 939 MeV is the average 
nucleon mass. V is the actual volume of the system, and 
Vf is so called free volume, which accounts for the finite 
size of nuclear species. We assume that all nuclei have 
normal nuclear density po, so that the proper volume of 
a nucleus with mass A is A/ p^. At low densities the 
finite-size corrections can be included via the excluded 
volume approximation, Vf/V ~ (l — p/po). This ap- 
proximation is commonly accepted in statistical models, 
and it is considered as a reasonable one at the densities 
p < O.lpo- Some information about the free volume at 
higher densities can be extracted from analysis of exper- 
imental data obtain ed in multifragmentation reactions 
(jScharenberg et al.1 (|2001D ). 

The internal excitations of nuclei play an important 
role in regulating their abundance, since they increase 
significant ly their entropy. Some authors (see, e.g., 
Rcf. (Ishiz uka. Ohnishi. fc Sumivoshi I (|2003[ ))) limit the 
excitation spectrum by particle-stable levels known for 
low excited nuclei. Within the SMM we follow quite 
different philosophy motivated by intensive investigation 
of disintegration of nuclei at temperatures more than 1 
MeV. Namely, we calculate internal excitation of nuclei 
by assuming that they have the same internal tempera- 
ture as the surrounding medium. In this case not only 
particle-stable states but also particle-unstable states 
will contribute to the excitation energy and entropy. This 
assumption can be justified by the dynamical equilib- 
rium of nuclei in the hot environment, and is supported 
by numerous comparisons with experiment (see section 
1). Moreover, in the supernova environment both the 
excited states and the binding energies of nuclei will be 
strongly affected by the surrounding matter. By this rea- 
son, we find it more appropriate to use an approach which 
can easily be generalized to include in-medium modifi- 
cations. The internal free energy of species {A, Z) with 
A > A is parameterized in the spirit of the liquid drop 
model, which h as been prov e d to be very successfu l in 
nuclear physics ()Bohr I (|1936D : IBondorf et ah I (|1995D ): 



Faz{T, p) = F^z + Fiz + FJ^ + F^z 



(5) 



Here the right hand side contains, respectively, the bulk, 
the surface, the symmetry and the Coulomb terms. 
The first three terms are taken in the standard form ( 



IBondorfet aTl (flggsh ). 
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where wq = 16 MeV, eq = 16 MeV, /3o = 18 
MeV, Tc = 18 MeV and 7 = 25 MeV are the 
model parameters which are extracted from nuclear phe- 
nomenology and provide a good description of multifrag- 
mentation data (Bondorf ct al. (1995); Botvina et ahl 

(1995; ); jScharenberg ct al (2 001); D'Agostin o et al. I 

([19961): iBellaize et al. I (l200l : lAvdevev et alTI ([2002); 
iD'Agostino et al. I (|1999D V However, these parame- 
ters, especially the symmetry coefficient 7, can be 
different in hot nuclei at multifragmentation condi- 
tions, and they should be determined from corre- 
sponding exper imental data (sec discussion in Refs. 
(IBotvina et al. I (12006): Buvukcizmcci et al. I pOOS) : 
iBotvina. Lozhkin. fc Trautmann I ()2002D )). 

In the electrically-neutral environment the nuclear 
Coulomb term should be modified to include the 
screening effect of electrons. This can be done, 
e.g., within the Wigner-Seitz approximation as wa s 
discussed in Re f . (iLamb. Lattimer. fc Pethick I (|1981); 
iLattimer et al. I ()1985[ )). We divide the whole system 
into spherical cells each containing one nucleus. The ra- 
dius of the cell is determined by the condition that it 
contains the same number of electrons as the number of 
protons in the nucleus. Then assuming a constant elec- 
tron density we obtain 



c(p) 
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2 VPop/ 



3 {eZf 



(9) 
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Pop J 



where rg = 1.17 fm and pop = {Z/A)po is the pro- 
ton density inside the nuclei. The screening function 
c{p) is 1 at pe = and at pe = Pop- To simplify 
calculations we oft en use an approx imation Pe/pov = 
p/pa, as in Ref. ([Lamb. Lattimer. fc Pethick (1981); 
ILattimer et al. I (|1985f )). which works well when neutrons 
mostly bound in nuclei, and leads to very similar results 
in many cases. We want to stress that the reduction 
of the Coulomb energy due to the electron screening fa- 
vors the formation of heavy nuclei. Nucleons and light 
nuclei (A < 4) are considered as structure-less particles 
characterized on l y by exact masses and proper volumes 
(jBondorf et al. I (J1995D ). Their Coulomb interaction is 
taken into account within the same Wigner-Seitz approx- 
imation. 

In Fig. [21 we present the pure nuclear pressure Pnuc 
(Eq. ([H)) caused by single nucleons and nuclei pro- 
duced after clusterization of nucleons. In addition, in 
the supernova environment free electrons existing to- 
gether with nuclear clusters modify the Coulomb en- 
ergy and lead to a negative contribution to pressure 
(|Bavm. Pethick. fc Sutherland ■ (|1971.) ). The physical 
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Fig. 2. — Comparison of results for the Coulomb contribu- 
tion Pc (bottom panels), pure Pnuc and total nuclear pressure 
Ptotnuc (top panels) as a function of temperature for Ye = 0.2, 
0.4 and p/po = 10"^ - 10"^ 

reason is that electrons compensate the positive charge 
of nuclei, and the Coulomb energy of Wigner-Seitz cells 
decrease when the density increases. This Coulomb con- 
tribution can be calculated as: 



AZ '^ 



(10) 



where the expression ([9]) for F^^ is used. The total nu- 
clear pressure Ptotnuc , which is the sum of the pure nu- 
clear pressure Pnuc and Pc, is demonstrated in Fig. [5] as 
a function of temperature. It is important that at T ^ 5 
MeV, when matter nearly completely dissociates into nu- 
cleons and lightest clusters, Pc is close to zero. In this 
region the total nuclear pressure coincides with the pure 
nuclear pressure. The Coulomb pressure becomes very 
important when heavy clusters dominate in the system. 
One can see in Fig. [5] that at low temperatures and high 
density the total nuclear pressure may be negative. In 
this case the nuclear clusterization is favorable for the 
collapse. However, the positive pressure of the relativis- 
tic degenerate electron Fermi-gas is considerably (more 
than an order of magnitude) larger. Therefore the to- 
tal pressure Ptot given by the sum of the total nuclear, 
electron and photon pressures in such environment will 
always be positive and the condition of thermodynamical 
stability ( dP*°^/dp > 0) will be fulfilled. 

As follows from Eq. Q, the fate of heavy nuclei de- 
pends strongly on the relationship between Faz and 
liAZ- In order to avoid an exponentially divergent con- 
tribution to the baryon density, at least in the thermody- 
namic limit {A — >■ oo), inequality Faz ^ Maz niust hold. 
The equality sign here corresponds to the situation when 
a big (infinite) nuclear fragment coexists with the gas 



of sm aller clusters (jBugaev. Gorenstein. fc MishustiiTI 
(|2001f )). When Faz > IJ'AZ, only small clusters with 
nearly exponential mass spectrum are present. However, 
there exists a region of thermodynamic quantities corre- 
sponding to Faz ~ fJ-AZ when the mass distribution of 
nuclear species is close to a power-law A^'^ with r w 2. 
This is a characteristic feature of the liquid-gas phase 
transition. The advantage of our approach is that we 
consider all the fragments present in this transition re- 
gion, and, therefore, can study this phase transition in 
all details. 

2.3. Leptons and photons 

Besides nuclear species the supernova matter contains 
leptons (e~, e"*", z/g, i^ei ■■■) and photons. At T, /i > nie 
the pressure of the relativistic electron-positron gas can 
be written as 
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where first-order correction (~ m^) due to the finite elec- 
tron mass is included, 5e=2 is the spin degeneracy fac- 
tor for electrons. The net number density pe and en- 
tropy density Se can be obtained from standard thermo- 



dynamic relations, 
which give 
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Electron neutrinos and anti-neutrinos are taken into ac- 
count in the same way, but as massless fermions, with the 
degeneracy factor twice smaller than for the electrons, 
i.e., gv=l. The photons are always close to the thermal 
equilibrium, and they are treated as massless Bose gas 
with zero chemical potential. The corresponding density 
M^, energy density e^, pressure P-y, and entropy density 



of photons are given by standard formulae: 



where g^—2. 
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3. THE SMSM EOS TABLES 

We have constructed the SMSM EOS ta- 
bles which are available on the Web at 

http://fias. uni-frankfurt. de/physics/mishus/research/smsm/ 



Below we will give the detailed information about mess 
of tables, entries and application of EOS tables. 

3.1. Mesh of the SMSM EOS Tables 
The SMSM EOS tables cover the following range: 

• Temperature: T = 0.2-25 MeV; for 35 T values. 

• Electron fraction Ye'. 0.02 — 0.56; linear mesh of 
Ye = 0.02, giving 28 Ye values. It is equal to the 
total proton fraction Xp, due to charge neutrality. 
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TABLE 1 

Temperature, electron fraction and density fraction 
values of smsm eos tables. 



TABLE 2 
SMSM EOS Tables are available in a machine-readable 

FORM IN THE ONLINE JOURNAL. 



T(MoV) 


0.20 


0.40 


0.60 


0.80 


1.00 


1.25 


1.50 


1.75 


2.00 


2.25 


2.50 


2.75 


3.00 


3.25 


3.50 


3.75 


4.00 


4.25 


4.50 


4.75 


5.00 


5.25 


5.50 


5.75 


6.00 


6.50 


7.00 


7.50 


8.00 


9.00 


10.00 


12.00 


15.00 


20.00 


25.00 


^e 


0.02 


0.04 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


0.22 


0.24 


0.26 


0.28 


0.30 


0.32 


0.34 


0.36 


0.38 


0.40 


0.42 


0.44 


0.46 


0.48 


0.50 


0.52 


0.54 


0.56 


P/PO 



lE-8 


1.78E-S 


3.16E-8 


5.62E-8 


lE-7 


1.78E-7 


3.16E-7 


5.62E-7 


lE-6 


1.78E-6 


3.16E-6 


5.62E-6 


lE-5 


1.78E-5 


3.16E-5 


5.62E-5 


lE-4 


1.78E-4 


3.16E-4 


5.62E-4 


lE-3 


1.78E-3 


3.16E-3 


5.62E-3 
lE-1 


lE-2 
1.78E-1 


1.78E-2 
3.16E-1 


3.16E-2 


5.62E-2 



• Baryon number density fraction p/ po = 

(10^^ — 0.32), giving 31 p/po values. 

The different T (35 points) , p/po (31 points) and Y^ 
(28 points) values sum up to 30380 EOS grid points are 
presented in the Table 1. We believe that the physics of 
SMSM EOS is sufficiently well presented with our grid 
points. Running time for calculations was approximately 
7.8 days (Work station:4xDual core AMD 2.2 GHz pro- 
cessor). We have obtained SMSM-Tables.txt, the size of 
this file 6.3 MB, it contains 31493 lines. 

3.2. Description of the SMSM EOS Tables 

The information is stored in a format which is 
very similar to the tables of IShen et al. I (119981) or 
IHempel fc Schaffner-BieliclTI (J2010D or iFurusawa et al. I 
( 20111 ) so that it can easilv be implemented in run- 
ning codes. The SMSM EOS tables file under the name 
"SMSM-Tables.txt" is the main file which is available in 
the online journal. 

In the beginning of file we write the general parameters 
of hquid-drop description: Wq = 16 MeV, Bq = 18 MeV, 
Tc = 18 MeV, 7 = 25 MeV, eo = 16 MeV, ro = 1.17 fm. 
The tables are written in the order of increasing T. In the 
tables, the values of T (MeV) are given at the beginning 
of each block, and the blocks with different T are divided 
by the string of characters cccccc. For each T, we present 
the results in the order of increasing Ye and p/po (see 
Table 1). We put also space lines between Yg blocks. 
A portion of the quantities in one line of SMSM EOS 
Tables are shown in Table 2 for guidance regarding its 
form and content. The quantities in one line of table for 
each density grid point the 19 different thermodynamic 
quantities are listed and defined as follows: 

1. Ye, electron fraction: Y^ — pe/ p , 

2. p/po, density fraction , 

3. Xn, neutron fraction: X„ = Pn/ P , 

4. Xp, proton fraction: Xp — Pp/p , 

5. Xd, doterium fraction: Xd = Pd/ P , 

6. Xt, tritium fraction: Xt = pt/ p , 

7. Xzfje^ ^He fraction: Xi^j^ — P^Hel P i 



Columnl 


Column2 


Column3 


Column4 


Column5 


Ye 


P/PO 


Xn 


Xp 


Xrf 


0.02 


O.lOOE-07 


0.9447 


0.3514E-38 


0.1248E-35 


0.02 


0.178E-07 


0.9445 


0.7178E-39 


0.4502E-36 


0.02 


0.316E-07 


0.9442 


0.1388E-39 


0.1609E-36 


Column6 


Column7 


Columns 


Column9 


ColumnlO 


^t 


^^H. 


^al-pha 


^h.a„„ 


< A > 


0.2067E~26 


a 


0.2873E-24 


0.05529 


89.66 


0.1327E-26 





0.6668E-25 


0.05553 


90.70 


0.8426E-27 





0.1519E-25 


0.05582 


91.76 


Columnll 


Columiil2 


Columnl3 


Columnl4 


Columnl5 


Etot 


Enua 


Stot 


Sntic 


ftot 


(MoVJ 


(McV) 


(ks) 


(kB) 


(McVfm~3) 


0.08460 


-0.1660 


13.81 


12.12 


0.4077E-09 


-0.02699 


-0.1666 


12.53 


11.58 


0.6277E-09 


-0.09008 


-0.1674 


11.57 


11.03 


0.1019E-08 


Columnie 


Column 17 


Columnl8 


Columnl9 




i'totr.uc 


Pe 


I'B 


f-g 




(MeVfm-3) 


(MoV) 


(MoV) 


(MoV) 




0.2S16E-09 


0.01716 


-2.063 


-16.96 




0.5011E-09 


0.03048 


-1.948 


-17.28 




0.8906E-09 


0.05394 


-1.833 


-17.60 





8. X alpha, alpha fraction X alpha = P alpha I P , 

9. Xheavy, heavy fragments {A > 4) frac- 
tion, -^ heavy Pheavyf P i 

10. < A >, average mass number of heaviest fragment , 

11. Etot, total energy is the sum of nuclear Enuc, 
electron E^ and photon Ephoton entropy contribu- 
tions. Etot = Enuc + Ee + En+ Ephoton, (McV/u) , 

12. Enuc, nuclear energy, (MeV/n) , 

13. Stot, total entropy is the sum of nuclear Snuc, 
electron Se and photon Sphoton entropy contribu- 
tions. Stot = Snuc + ie + Sn + Sphoton, (^ij) , 

14. Snuc, nuclear entropy, (fcs) , 

15. Ptot, total pressure is the sum of the pressure of 
nuclear species -Ptotnuc, electrons Pc, neutrons P^ 
and photons Pphoton, (MeVfm^^) , 



Pt, 



Pu 



I ^c I -fii I .fphoton , 

16. Ptotnuc, total nuclear pressure, (MeVfm^'^) , 

-M:otnuc — .'nuc i .t C i 

17. Pe, chemical potential of electrons, (MeV) , 

18. p,B, chemical potential of mass number, (MeV) , 

19. pq, chemical potential of charge, (MeV) . 

3.3. Application of the EOS tables and illustrative 
results 

In addition to standard output (the 19 different ther- 
modynamic quantities), like the pressure, entropy energy, 
fr actions of iew light p articles and heavy nuclei (actually 
in iShen et al. I ()1998D case a single heavy nucleus) , we 
provide a possibility to extract the ensemble of all heavy 
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nuclei. We are doing it by providing the chemical poten- 
tials and the formula 






3 



(15) 



exp 



-7^{FaZ- f^Az) 



which is general for all points, and which can calculate 
densities of nuclei with all mass numbers and charges. 
We believe that in simulations it could be quite helpful 
to have knowledge about all heavy nuclei at each points, 
since weak reactions (with neutrinos and electrons) can 
be calculated much more precisely with the ensemble of 
heavy nuclei. Thus, one can easily calculate the the en- 
semble of all heavy nuclei of stellar matter using this 
formula and table parameters. 

The mass distributions contain important informa- 
tion about fragmentation of matter in the nuclear 
liquid-gas phase transition (coexistence region). The 
concept of statistical equilibrium assumes a continu- 
ous interaction between fragments via specific micro- 
scopic processes, like absorption and emission of neu- 
trons, which provide equilibration ( see, e.g., discussion 
in Ref. (jBotvina fc MishustiiTI ()201Cir ))). In this situation 
the nuclei can remain hot and have modified properties 
and masses, which are different from the cold isolated 
nuclei. We should note that it is stored also "SMSM- 
Relative-Mass-Yield.dat" data file for mass distributions 
of system for each calculations to make quick check one 
by one for us. This file size is 356 MB, it contains 30380 
lines and every line has 1000 column (A=l-1000). Since 
we have 35 T and 28 Y^ values, (35(T) x 28 (ye)==980 val- 
ues for 31 p/po values), we grouped 4 density intervals 
from lower to higher as 



Density- l=/3//7o = 10 



.62.10- 



Density-2=p/po = 10"*^ - 5.62.10" 



values) 



values) , 



• Density-3=p//9o = lO"'^ - 5.62.10"^ (8 values), 

• Dcnsity-4=p/po = 10"^ - 3.17.10"^ (7 values). 

Using these interval for data, we have plotted 3920 
figures as eps files which contains information about 
mass distributions of system for each grid point. We 
have presented example figures for T = 1, and 3 MeV, 
Ye = 0.2 - 0.4 and different densities in Figs. H El [5] and 
[HI As seen in these figures, we also put Xn, Xp, Xaipha, 
Yheavy, < A >, pa{= Pb) and pz{^ Pq) values from 
tables with the same order of p/po as the control param- 
eters inside of figures. These figures are presented on the 
SMSM web page with 

• SMSM-Massdist-Figs-Density-l.zip, 

• SMSM-Massdist-Figs-Density-2.zip, 

• SMSM-Massdist-Figs-Density-3.zip, 

• SMSM-Massdist-Figs-Density-4.zip 

files, each file contains 980 eps figures. For example, in- 
side SMSM-Massdist-Figs-Density-4.zip file, one can find 
mass distribution figure for T = 1 MeV, Ye = 0.20 and 
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Fig. 3. — Mass distributions of fragments produced in matter 
with temperatures T = 1 MeV, electron fractions Ye — 0.2 
and densities p/pa = 10~^ — 3.17.10"^. (Color version online.) 
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Fig. 4.— The same as Fig.[3]but at Y^ = 0.4. 



'f.eps 



p/po = 10-2-3.17.10-1 in T-l-Ye-020-DENSITY-. 
file as shown in Fig. |31 

At very low temperatures the SMSM predicts a 
Gaussian-like distribution for heavy nuclei as seen in 
Figs. [3] and ID In this case an approximation of a 
single heav y nucleus ado pted in the Lattimer fc Swestyl 
(|1991h and iShen et al. I (1998. ) EOSs may work reason- 
ably well for calculations of thermodynamical character- 
istics of matter. However, already at T > 1 MeV the gap 
between the Gaussian peak and light clusters and nucle- 
ons is essentially filled by nuclei of intermediate masses 
leading to characteristic U-shaped distributions. The ex- 
amples of such distributions are shown in Figs. [5]and|6l 
They give a typical evolution of mass distributions of nu- 
clei in the coexistence region of the liquid gas phase tran- 
sition. Previously it wa s demon s trated for disintegration 
of finite nuclei (.Bondorf et al. I (|1995l) ). These distribu- 
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Fig. 5.— The same as Fig. |3] but at T = 3 MeV. 
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Fig. 6.— Tiie same as Fig. [5] but at Ye = 0.4. 

tions continuously transform with decreasing density as 
well as with increasing temperatures into exponential fall 
of fragment yields with A. 

Information about isotopic yields helps to understand 
the differences observed for summed quantities like mass 
yields or mass fractions. It is also possible to obtain 
isotopic distributions of stellar matter using SMSM EOS 
tables and Eq. (J15p . Here, we have shown examples from 
isotopic distributions oi Z = 8 and Z = 26 in Fig. [71 As 
seen from these figures, the SMSM gives Gaussian type 
distributions, which are the consequence of the liquid- 
drop description of fragments. 

Recently, we have compare d mass and isotope 
distributions of SMSM (iBotvina fc MishustiiTI 



(l20TnhV iHemne fc Schaffner-Bielich I (120101 1 and 



ipe & bcJ 

iFurusawa et al. I ()2011[ ) EOS models at subnuclear 
densities. Ye = 0.2,0,4, temperatures (0.5-10 MeV) 
expected to occur during the coll apse of massive stars 
and supernova explosions in Ref. iBuvukcizmcci ct alTI 




50 60 70 80 90 100 
A 



Fig. 7. — Isotopic distributions of Z = 8 and Z — 26 frag- 
ments produced in matter with temperatures T = 1,2,3, 
and 5 MeV, electron fractions Ys = 0.2 and 0.4, densities 



p/po = 10- 



10" 



(j2013f ). presently, we have found significant differences 
between mass distributions of SMSM, HS and FYSS 
models, which use different assumptions on properties 
of fragments, especially at low electron fractions, low 
temperatures and high densities. We note that these 
differences are also shown in the behaviour of nuclea r 
pressure, see Fig. 22 in Ref. iBuvukcizmeci et al. I (|2013| ). 

3.4. Nuclear composition under condition of 
l3- equilibrium in stellar matter 

An important case when full /3-equilibrium can be 
reached in matter (see Eq. ^). One can expect it in 
some processes like slow collapse of massive stars, in 
crust of neutron stars, and in some conditions of stel- 
lar matter after supernova explosions. On the other 
hand, the /^-equilibrium can be a useful physical limit for 
theoretical estimate of the nuclear composition without 
full knowledge of weak rea ctions. The SMSM all o ws for 
this kind of calculations (iBotvi na fc Mishustin I ()2004 
l2010f )) and the corresponding Tables, at baryon densi- 
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Fig. 8. — Electron fractions in /^-equilibrium case versus tem- 
perature at densities p/po = 0.3 — 10~^. 



ties p = (10^^ - 0.32)po, and temperatures T = 0.2 - 25, 
are under construction. In this case the corresponding 
electron fraction Ye can be systematically found at given 
p and T. In the following this Ye may be used in stan- 
dard SMSM tables too. In Fig. |8]we show the calculated 
electron fraction under /3-equilibrium condition for some 
densities versus temperature. One can see an interesting 
behavior: There is a trend of increasing Yg at large tem- 
peratures, which has a simple explanation: At high T the 
constraint imposed on electrons by their large chemical 
potentials (at high densities) becomes weaker. At the 
largest density (p = 0.3po)j when heavy nuclei always 
present in the ensemble at all temperatures, this trend is 
permanent. However, at smaller densities (p ^ O.Olpo) 
the heavy nuclei can exist only at low temperatures since 
they disintegrate into nucleons and light clusters at high 
temperatures (T > 3 — 5 MeV). These nuclei have big 
charges and present in a large amount, that requires also 
a large number of electrons at large temperature. As a 
result, at such densities there is a temperature at which 
the electron fraction has a minimum, corresponding to 
transition from heavy nuclei to free nucleons and light 
clusters. This is a very instructive effect of coupling of 
nuclear composition and electron fraction in matter. 



4. CONCLUSIONS 

We have presented the constructed SMSM EOS tables 
which give comprehensive information on physical prop- 
erties of stellar matter at subnuclear densities, temper- 
atures T < 25 MeV and electron fractions between 0.02 
and 0.56. The SMSM EOS can be used for hydrodynam- 
ical simulations of massive stars collapse and supernova 
explosion processes. We describe the whole ensemble of 
nuclei, without any artificial limitation of their masses 
and charges, and general characteristics of matter, like 
temperature and chemical potentials. The nuclear mass 
and isotope distributions are needed for realistic calcula- 
tions of electron capture and neutrino induced reactions. 
As we have noted the present EOS at large densities 
as directly related to nuclear multifragmentation reac- 
tions, and the properties of nuclei embedded in dense 
surrounding can be taken by obtained from analysis of 
these reactions. Our calculations of the thermodynami- 
cal parameters and nuclear composition of matter is espe- 
cially instructive in the case of /3-equilibrium. A special 
case of /3-equilibrated nuclear matter is considered too. 
In particular, we have found non-monotonous behavior 
of the electron fraction as a function of temperature at 
p < O.lpo- This effect is connected with the rearrange- 
ment of the fragment mass distribution from U-shape to 
exponential fall-off. These results can be used for model- 
ing proto-neutron stars and crusts of neutron stars. We 
believe that using EOS developed by different groups of 
authors is crucially important for evaluating real behav- 
ior of stellar matter. 

5. ACKNOWLEDGEMENT 

N.B. and A.S.B. are grateful to Frankfurt Institute 
for Advanced Studies (FIAS) for support and hospi- 
tality. This work is supported by the Helmholtz In- 
ternational Center for FAIR within the LOEWE pro- 
gram. N.B. thanks Selcuk University-Scientific Research 
Projects (BAP) for partly support. N.B. acknowledges 
that a part of the numerical calculations were carried 
out on the computers at Physics Department in Selcuk 
University. N.B. thanks to R.Ogul for encouragement, 
A.E.Kavruk and T.Ozturk for the help for the usage of 
computers in S.U.. A.S.B. acknowledges the support by 
the Research Infrastructure Integrating Activity 'Study 
of Strongly Integrating Matter' HadronPhysicsS under 
the 7th Framework Programme of EU. I.N.M. acknowl- 
edges support from the grant NSH-2015.2012.2 (Russia). 



REFERENCES 



Arcones, A., Martmez-Pinedo, G., O'Connor, E. et al. 2008, 

Phys. Rev. C, 78, 015806 
Avdeyev, S. P., Karnaukhov, Y. A., Petrov, L. A. et al. 2002, 

Nud. Phys. A, 709, 392 
Baym, C, Pethick, C, & Sutherland, P. 1971, ApJ, 170, 299 
Bellaize, N., Lopez, C, "Wieleczko, J. P. et al. 2002, Nud. Phys. 

A, 709, 367 
Bohr, N. 1936, Nature, 137, 344 

Bondorf, J. P., Botvina, A. S., Iljinov, A. S., Mishustin, I. N., & 
1995, Phys. Rep., 257, 133 



Sneppen, K 
Botvina, A. S., Iljinov, A. S. 

Nud. Phys., 42, 712 
Botvina, A. S., Iljinov, A. S. 

Phys. A, 507, 649 
Botvina, A. S., Mishustin, I. N., 

1995, Nud. Phys. A, 584, 737 
Botvina, A. S., Lozhkin, O. Y., & Trautmann, W, 

Rev. C, 65, 044610 
Botvina, A. S. & Mishustin, I. N. 2004, Phys. Lett. B, 584, 233 



& Mishustin, I. N. 1985, Sov. J. 
& Mishustin, I. N. 1990, Nud. 
Begemann-Blaich, M. et al. 

2002, Phys. 



Erdogan, M. et al. 2006, 



2010, Nud. Phys. A, 843, 98 
& Mishustin, I. N. 2001, Phys. 

Mishustin I. N., & Ogul, R. 



Botvina, A. S., Buyukcizmeci, N. 

Phys. Rev. C, 74, 044609 
Botvina, A. S. & Mishustin, I. N 
Bugaev, K. A., Gorenstein, M. I. 

Lett. B, 498, 144 
Buyukcizmeci, N., Botvina A. S.. 

2008, Phys. Rev. C, 77, 034608 
Buyukcizmeci, N., Botvina, A. S., Mishustin I. N. et al. 2013, 

Nud. Phys. A, 907, 13 
D'Agostino, M., Botvina, A. S., Milazzo, P. M. et al. 1996, 

Phys. Lett. B, 371, 175 
D'Agostino, M., Botvina, A. S., Bruno, M. et al. 1999, Nud. 

Phys. A, 650, 329 
Dynamics and Thermodynamics with Nuclear Degrees of 

Freedom, ed. by Chomaz, Ph., Gulminelli, F., 

Trautmann, W., & Yennello, S. J. (Springer, 

Berlin/Heidelberg/New York) 2006, Eur. Phys. J. A, 30 



Tabulated equation of state for supernova matter 



Furusawa, S., Yamada, S., Surniyoshi, K., & Suzuki, H. 2011, 

ApJ, 738, 178 
Gross, D. H. E. 1990, Rep. Prog. Phys., 53, 605 
Hauger, J. A., Srivastava, B. K., Albergo, S. et al. 2000, Phys. 

Rev. C, 62, 024616 
Hempel, M. & Schaffner-Bielich, J. 2010, Nud. Phys. A, 837, 210 
Hempel, M., Schaflfner-Bielich, J., Typcl, S., & Ropkc, G. 2011, 

Phys. Rev. C, 84, 055804 
Hempel, M., Fisher, T., Schaffner-Bielich, J., & Liebendorfer, M. 

2012, Astrophys. J., 70, 748 
Hudan, S., Mcintosh, A. B., Black, J. et al. 2009, Phys. Rev. C, 

80, 064611 
Iglio, J., Shetty, D. V., Ycnnello, S. J. et al. 2006, Phys. Rev. C, 

74, 024605 
Ishizuka, C., Ohnishi, A., & Surniyoshi, K. 2003, Nucl. Phys. A, 

723, 517 
Lamb, D. Q., Lattimcr, J. M., Pcthick, C. J., & Ravcnhall, D. G. 

1978, Phys. Rev. Lett., 41, 1623 
Lamb, D. Q., Lattimer, J. M., & Pethick, C. J. 1981, Nucl. 

Phys. A, 360, 459 
Lattimer, J. M., Pethick, C. J., Ravcnhall, D. G., & 

Lamb, D. Q. 1985, Nucl. Phys. A, 432, 646 
Lattimcr, J. M. & Swcsty, D. F. 1991, Nucl. Phys. A, 535, 331 



Ogul, R., Botvina, A. S., Atav U. et al. 2011, Phys. Rev. C, 83, 

024608 
Pienkowski, L., Kwiatkowski, K., Lefort, T. et al. 2002, 

Phys. Rev. C, 65, 064606 
Prakash, M., Bombaci, I., Prakash, M. et al. 1997, Phys. Rep., 

280, 1 
Rodionov, V., Avdcyev, S. P., Karnaukhov, V. A. ct al. 2002, 

Nucl. Phys. A, 700, 457 
Scharenberg, R. P., Srivastava, B. K., Albergo, S. ct al. 2001, 

Phys. Rev. C, 64, 054602 
Shen, H., Toki, H., Oyamatsu, K., & Surniyoshi, K. 1998, Nucl. 

Phys. A, 637, 435 
Souliotis, G. A., Botvina, A. S., Shetty, D. V. ct al. 2007, 

Phys. Rev. C, 75, 011601 
Sugahara, Y. & Toki, H. 1994, Nucl. Phys. A, 637, 557 
Sumiyoshi, K., Yamada, S. Suzuki, H. ct al. 2005, ApJ, 629, 922 
Surniyoshi, K. & Ropkc, G. 2008, Phys. Rev. C, 77, 055804 
Typel, S., Ropke, G., Klahn, T., Blaschkc, D., & Woltcr, H. H. 

2010, Phys. Rev. C, 81, 015803 
Viola, v., Lefort, T., Beaulicu, L. et al. 2001, Nucl. Phys. A, 

681, 267 
Wang, G., Kwiatkowski, K., Bracken, D. S. et al. 1999, 

Phys. Rev. C, 60, 014603 



